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Abstract 

We formulate a dynamical model for microtubules interacting with a catastrophe-inducing 
boundary. In this model microtubules are either waiting to be nucleated, actively growing or 
shrinking, or stalled at the boundary. We first determine the steady-state occupation of these 
various states and the resultant length distribution. Next, we formulate the problem of the Mean 
First-Passage Time to reach the boundary in terms of an appropriate set of splitting probabilities 
and conditional Mean First-Passage Times, and solve explicitly for these quantities using a differen- 
tial equation approach. As an application, we revisit a recently proposed search-and-capture model 
for the interaction between microtubules and target chromosomes [Gopalakrishnan &: Govindan, 
Bull. Math. Biol. 73:2483-506 (2011)]. We show how our approach leads to a direct and compact 
solution of this problem. 
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I. INTRODUCTION 



Microtubules (hereafter abbreviated as MTs) are filamentous macro-polymers built from 
tubulin dimers. They are one of the components of the cytoskeleton of all eukaryotic cells. 
They play a number of roles ranging from providing mechanical stability to the cell, serving 
as transport pathways enabling linear transport of versicular cargo by motor proteins, and 
providing forces for the positioning of organelles and other intracellular components (for a 
general overview see pj]). Perhaps their most striking function appears during mitosis, where 
they form the mitotic spindle, the machinery for positioning and separating the duplicated 
chromosomes prior to cell division. They owe their functional plasticity to an intrinsic so- 
called 'dynamical instability' mechanism j2| that causes individual MTs to stochastically 
alternate between growing and shrinking states. By controlling this dynamical process, 
through MT associated proteins (MAPS) that nucleate new MTs or selectively stabilize or 
destabilize them by locally or globally changing the rates with which they switch between 
dynamical states, cells are able to reconfigure MT assemblies on timescales as fast as a 



few minutes. The canonical mode^ 
1990's by Dogterom and Leibler 



to describe MT djTiamics was developed in the early 
3|. This model showed that isolated MTs, depending on 
their dynamical parameters, can either be in a regime of bounded growth leading to an 
exponential length distribution in a steady state, or in a regime of unbounded growth in 
which the MT length on average increases linearly in time. Of course, MTs 'live' within the 
confines of a finite size cell, whose dimensions (~ 10/im) are comparable to the observed 
lengths of MTs. Interactions between MTs and boundaries, be it the cell cortex or the 
surface of other intracellular compartments, are therefore important. Indeed, a number of 
MT functions depends critically on these interactions: examples are nuclear positioning in 

yeast [4], spindle positioning in C. elegans [5;], and the orientation of the cortical MT array 

I 

in plant cells |6j. 

In spite of this clear relevance, it appears that a systematic approach to the theory of 
MTs interacting with boundaries is lacking from the literature. The one problem of this 
type which did receive substantial attention, is the search-and-capture mechanism by which 
MTs are thought to find the condensed chromosomes prior to mitosis [rljsljQ] 10 1, which 
involves estimating the mean first-passage time of a MT to hit a limited size target at a 
distance from its nucleation point. However, although these works in fact do contain some 
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of the basic features of the MT-boundary problem, it is mostly hidden (literally in the case 
of Ref. js], actually an unpublished thesis) under the specifics of the intended application. 
Moreover, these works also rely heavily on 'forward' techniques involving the time evolution 
of the full probability density for reaching a given state from specified initial conditions. 
Although this approach obviously yields a full solution of the problem, for passage-time 
problems, which effectively require integrating over final states, the full probability density 
is in a sense a form of "overkill". The treatment of this type of problems can in fact be 
simplified considerably by using 'backward' techniques, as is e.g. elegantly illustrated for 



diffusion problems in Redner's monograph 



In the present work we show how this 



approach can be used from the ground up to solve the problem of a MT interacting with a 
boundary. 

The outline of the paper is as follows: In Section [III we introduce the dynamical model of 
a MT interacting with a boundary, solve for its steady state properties, such as the average 
length, and choose an appropriate set of non-dimensionalized parameters and variables. 
In Section IIIII we turn to the analysis of the mean first-passage time to the boundary, 
formally solving this in terms of a small set of splitting probabilities and conditioned mean 
first-passages times, which are subsequently determined explicitly through the solution of 
appropriate linear boundary value problems. We then use some biological data on MT 
dynamics to estimate order of magnitudes for the quantities involved. Finally, in Section IIVI 
we revisit the recent search-and-capture model discussed by Gopalakrishnan and Govindan 
|lol |. and show how it is compactly solved using the techniques introduced. We then finish 
with a number of concluding remarks and two technical appendices. 



II. MEAN LENGTH 
A. Dynamical model 

The standard two-state dynamical instability model describes MTs with length / that are 
either growing with velocity f+ or shrinking with velocity f _ and can switch between the 
growing and the shrinking state (a catastrophe) with a constant rate r+ and between the 
shrinking and the growing state (a rescue) with a constant rate r_. Collectively we call these 
two states the active states, and denote the corresponding state space by A. We extend this 
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FIG. 1: Schematic of our model showing the three components of the state space: The nucleation 
state N from which MTs are nucleated at rate r„. The active state A, in which MTs either grow 
or shrink, with velocities and V-, and switch between these two states with rates r-|_ and r_, 
and a barrier state B from which MTs exit at a rate rf,. 



model by two more states: a nucleation state A^, in which a MT enters upon shrinking back 
to zero length, and from which it can be (re) nucleated into a zero- length growing state at a 
constant rate r„, and a boundary state B, in which a MT enters upon hitting a boundary 
at a distance L from the nucleation point, and which it leaves in a shrinking state at a rate 
rf,. Formally the state space of this extended model is therefore given by f2 = U A U 5. 
We illustrate the model and its state space in Figure [1] It should be noted that, from a 
biophysical point of view, the model is of course an idealization. In reality a growing MT 
impinging on a boundary will generate forces. These forces will affect the growth velocity 
and the propensity^ to switch to the shrinking state, so that the latter is no longer a simple 
Poisson process [12|]. Moreover, these forces may also deform the boundary, lead to buckling 



of the MT, or cause it to slide along the boundary (for a review see [iS]). These additional 
complexities, however, will have limited impact on the results to be presented here, as long 
as the residence time at the boundary is small compared to the time to traverse the distance 
between nucleation point and boundary, which is certainly the case for effectively reflective 
boundaries for which r;, 3> r+, and we choose to ignore them here. 

The dynamics of the model described above defines a time-homogeneous Markov process 
on the full state space. It is, however, technically convenient to split the dynamics on the 
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'active' part of the state space A from that on the two 'waiting' states N and B, and deal 
with the communication between these different states through boundary conditions. We 
therefore first define the probabihty densities (per unit length) (/, tjcuo, to) for an active 
MT to have length / and be in state s = +, — at time t, given that it was in some state cjq 
at time to <t. These densities satisfy the evolution equations 

d d 

— m+ (/, t\LJo, to) = -v+—m+ (/, t\LJo, to) - r+m+ (/, t\LJo, to) + r_m_ (/, t\uo, to) , (1) 
ot ol 

d d 

— m_ {I, t\uJo, to) = v-—m- {I, t\uJo, to) - r-m^ {I, t\uJo, to) + r+m+ {I, t\uJo, to) . (2) 

Likewise we define the probability M„ {t\uJo,to) for the MT to be in the nucleation state at 
time t, given that it was in some state a;o G at time to < t. This probability satisfies the 
evolution equation 

d 

-Q^Mn {t\uo, to) = v^m^ (0, t\uo, to) - r^M^ {t\uo, to) . (3) 

The probability Mb {t\uio-i to) for the MT to be in the boundary state at time t in turn satisfies 

d 

—Mb {t\uJo, to) = v+m+ {L, t\uJo, to) - nMb {t\uJo, to) (4) 
This system of equations is closed by the boundary conditions 

v+m+ (0, t\uJo, to) = r,,Mn (^l^o, ^o) , (5) 
v-m- (L, t\uJo, to) = nMb {t\uJo, h) ■ (6) 

By construction the dynamics on the full state space conserves probability, and indeed if we 
define the total probability 

M {t\uJo, to) = Mn (tjcUo, to) + Ma (t|cJo, to) + Mfe (tjcjo, to) 

poo 

= M„ (t|a;o, to) + / dl {m+ (Z, t|a;o, to) + m_ (Z, t|a;o, to)} + M^ (t|a;o, to) (7) 

then 

— M(t|a;o,to) = 0, (8) 

allowing us to set M (t\uJo,tQ) = 1. 
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B. Steady state behaviour 



In steady state the probabilities do not depend on time nor on initial conditions, allowing 
us to write the evolution equations as a set of balance equations 

f+— m+ (/) = —r^m^ (/) + r^m^ (/) , (9) 
—V-—m_ (/) = —r_m_ (/) + r_|_m_|_ (/) , (10) 

v_m_ (0) = r„M„, (11) 
v+m^ (L) = r^Mfe, (12) 

to be supplemented by the boundary conditions 

v+m+ (0) = r„M„, (13) 

v_ra_ (L) = TbMfe. (14) 

Adding Eqs. (ED and <^ yields 

|Km+(0-t;-m_(0} = 0. (15) 

Combining either Eqs. f lTT]) and f lT^ . or Eqs. f lT^ and dH]), shows that the constant of 
integration vanishes, and hence 

v+m_|_ {I) = t>_m_ (/) . (16) 

This identity is now used to eliminate m_ (/) from from which we then readily find that 
for / < L 

m+ (/) = !:!^e-'/'", (17) 
m_ (/) = !i^e-'/\ (18) 

where the length 

(19) 

is of course only positive when r^V- — r_f + > 0, the so called bounded-growth regime, and 
represents the steady-state average length of an active MT in the absence of the boundary. 
Although in the presence of boundaries one can also consider the unbounded growth regime 



r_|_f _ — r^v^ < 0, as was done e.g. in [l^ for the case where the boundary is fully reflecting, 
we will not do so here, and focus exclusively on the bounded growth case. The dependence 
of the probability that the MT is in the boundary state can, from fll2p . be shown to obey 
the following relationship 

Mb = —e-TMn (20) 
rb 

The normalization condition ([7]) can then be used to determine the probabilities for an MT 
to be in a nucleation state, an active state and a boundary state respectively 

M„ = ^ , (21) 

l + rj(^ + ^) (1 - e-Vr) + r.e-Vr 



l + rj(^ + ^) (1 - e-Vr) + r.e-vr 



Ma = ' --' (22) 



r-n Q-L/l 

Mb = y (23) 

l + rj(^ + ^) (1 - e-Vr) + r.e-Vr 



X 



Taking the limit L — )■ oo we find Ma oc r„/ + As density oc nucleation rate 

lifetime, this suggests that the time 

f.r(±,i).^±±^ (24) 

is the expected lifetime of a, otherwise unconstrained, length-zero newly nucleated MT, a 



result indeed first derived by Rubin 15 1. 



We now define the mean length of the active MTs 

f 1 _ (1 + I.) e-V'" 



(^)a = ^J^ dll{m^{l)+m.{l)} = l ^ (1-e-Vr) ' ^^^^ 

The time-averaged length of the MT over the full ensemble, is then simply 

(/) = Ma + M,L (26) 

We can readily check the limits (/)^ — )• / when L — i- oo, and (l)^ ^ |L for L — )■ 0. The 
latter limit can be understood by considering that for very small L the growing and shrinking 
traversal times L/v+ and L/v- become small with respect to the mean time between between 
catastrophes, l/r+, and rescues, l/r_, respectively, so that the MT is deterministically 
"bouncing" between the endpoints / = and I = L. 



Symbol 


Definition 


Meaning 




1!+ 

V 


growth speed 




V 


shrinkage speed 


P+ 


r^t 


catastrophe rate 


P~ 




rescue rate 


Pn 


rj 


nucleation rate 


Pb 


nt 


barrier unbinding rate 


A 


L 
I 


distance to barrier 



TABLE I: Dimensionless parameters of the MT dynamical model 
C. Dimensional analysis 



We have deliberately deferred the dimensional analysis of the system up to this point to 
allow the results of the steady state solution to guide us to a natural choice of the units of 
length and time. In view of fl2T|) and ( l25l) . we choose I ([21]) as the unit length, and i as the 
unit of time. For completeness sake, we can also introduce the unit of speed 

1 i 



1 1 



(27) 



By convention we will adopt the Greek alphabet do denote dimensionless quantities. We 
can now introduce the dimensionless parametersNote that this assignment, which has the 
clear advantage of maximizing the interpretability of the non-dimensional equations, does 
have the disadvantage of leaving dependencies among the parameters, as by construction 

1, (28) 



1 1 
— + — 



(29) 

To denote the independent variables of time and length we write r and A respectively. Finally, 
the densities, as our dependent variables, are denoted by /i^ (A, t\ujo, tq) = Inis [XI, riluo, rot) . 
Using these notations the steady-state results of the previous section can be summarized 

as 

1 



l + p„(l-e-^) + ^e-^' 



Pb 



(!-(! + A) e-^) 



(1 



(30) 
(31) 
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For future reference we will also rewrite the evolution equations in a more compact 
notation. To do so we treat pairs of functions (v5+(A, r), ip-{X, r)), defined on the growing and 
shrinking parts of the state space respectively, as a single vector valued function ipsi^, ''")) -s = 
+, — . This allows us to write 



The fact that we use the notation G*, signifying the Hermitean conjugate of the generator 
G of the Markov process, is conventional when discussing the forward Kolmogorov equation, 
which is the formal term for the evolution equation for the probability densities jl^. 

III. THE MEAN-FIRST PASSAGE TIME 

We approach the problem of calculating the Mean-First Passage Time (MFPT) for a 
microtubule to hit the boundary at a distance L in three steps. We first provide a formal 
solution to the problem in terms of suitably chosen set of survival (and ruin) probabilities. 
We then calculate the static splitting probabilities that describe the relative weights of the 
direct and indirect paths of reaching the boundary, and finally determine the conditional 
MFPTs corresponding to these sets of paths. 

A. Formal solution 

We first define the survival set, the subset of state space excluding the boundary state, 
i.e. = Q/B = N U A. Our goal is to determine the survival probability of the process 
in this set starting at r = from an arbitrary active state with length A < A, which 
we denote by Sq^ (r|A, s). We will deconstruct this survival probability with the aid of a 
number additional probabilities defined as follows: (r|A, s), the probability of surviving 
in the active part of the state space, i.e. of not having passed either A or shrunk back to 
zero length from the initial condition, and the conditional ruin probabilities (t|A, s) and 




(32) 



where G* is the operator matrix 




(33) 
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(r|A, s), being the probabilities to have exited into the boundary state at A = A or the 
nucleation state at A = at time r respectively, without leaving the active state at any 
prior moment. Finally, we define the survival probability to remain in the nucleation state, 
which, because nucleation is a simple Poisson process, is given hj ('^1^) = 6xp {—pnr). 
Each of these survival or ruin probabilities has a corresponding waiting time distribution, 
symbolically given by a {t\u) = —^5' {t\u) or a {t\u) = -^R {t\u). 
We now note the following two identities 

Sn^ (r|A, s) = Sa (r|A, s) + [ dr'a^ (r'|A, s) Sn^ (r - t'\N) (34) 

Jo 

Snj, {r\N) = Sjv (r|iV) + [ dr'a^ (r'|iV) S^^ (r - r'|0, +) (35) 



^0 

The first captures the fact that starting from an active state the MT survives either by 
remaining active and not reaching A, or by shrinking back to zero at some intermediate 
time and then surviving from the nucleation state. The second states that a MT in the 
nucleation state survives either by remaining in this state, or being nucleated into a growing 
one at an intermediate time and then surviving from the zero-length growing state. As 
survival in the active state means not exiting either at length A into the B state or at length 
into the N state we have that 

Sa (r|A, s) = 1 - i?f (r|A, s) - R^ (r|A, s) (36) 

As is clear from the steady state solution a MT will always leave the active state for large 
enough time (M^ > independent of the initial conditions) so that Sa (oo|A, s) = 0. The 
ultimate conditional ruin probabilities R^ {oo\X, s) and R^ (c>o|A, s) are usually, and aptly, 
called splitting probabilities, as the total ruin probability is 'split' betwen them R^ (oo|A, s) + 
R^ {oo\\,s) = 1. We can thus rewrite identity ( l34|) as 

Sn^ (r|A, s) = [R^ (oo|A, s) - R^ (r|A, s)) + « (oo|A, s) - R^ (r|A, s)) 

+ r dr'a^ (r'|A, s) Sn^ (r - r'\N) (37) 

We can now define 

/•oo /"OO 

Tn^{X,s)= / dr Snj^{T\X,s) = / rfrr (Tq^ (r|A, s) (38) 
Jo Jo 
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as the MFPT for the process to pass length A starting from the active state (A, s). Integrating 
the first two terms on the righthand side of (1371) over time yields 



roo rod 

/ dr {R^{oo\X,s) - R^{t\\s)) = dr 
Jo Jo 



oo 

Ta^ (r|A,s) 



= Rl (oo|A, .) \^ / \ ' / ^ Rl (oo| A, .) (A, .) (39) 
Jo dTa^[T\\s) 

which introduces the conditional MFPT of the process to exit at A without ever shrinking 
to 0, and similarly 

poo 

/ dr {rI{oo\\s)-RI{t\\s))=RI{oo\\s)T^{\s) (40) 
Jo 

Integrating the final term gives 

/ dr / dr'a'l (r'|A, s) Sn^ (r - r'\N) = 
Jo Jo 

poo poo 

/ dT'a^{T'\X,s) dTSnAr-r'\N) = RAioo\X,s)TnAN) 

Jo Jt' 

The MFPT from the nucleation state is readily obtained from (l35ll and yields 

Tn, (N) = (N) + Tn, (0, +) = ^ + Tn, (0, +) (41) 

Pn 

where we have used that exiting the nucleation state is sure i.e. Sn (oo|A^) = 0. Collecting 
all these results then yields 

Tn^ (A, s) = R^ (oo|A, s) Tf (A, s) + R^ (oo|A, s) (A, s) + - + Tn^ (0, +)1 (42) 



k Pn J 

The interpretation of this result is clear: Starting from (A, s) the MT either exits directly at 
A = A, which happens with probability R^ (oo|A, s) and (on average) takes time Tf (A, s), 
or the MT first shrinks back to A = (with probability i?;^ (oo|A, s)) which takes a time 
(A,s), and then has to wait a time to be renucleated, after which it takes time 
Tn^ (0, +) to reach A = A for the first time. We illustrate this result schematically in Figure |2j 

If we now consider a MT starting at A = and in the growing state s = + , we can 
selfconsistenly solve for the MFPT Tq^ (0, +), which is our main result 



(0. +) = (0, +) + ^11^ + ^ ^ («) 
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X 



A 



Splitting 




+ 



FIG. 2: Schematic showing how the unconditional MFPT (grey arrow) for a growing MT to reach 
the boundary from starting length A splits into two conditional MFPTs: The first associated with 
all direct paths (black arrows) from A to the boundary at A. The second associated with all 
paths that shrink to zero length without reaching the boundary, spend time in the nucleation state 
(wiiglly line) and then tries again unconditionally growing from length 0. 



This latter result is in fact sufficient to solve the general problem, because in order to reach 
the boundary at A from length the MT first has to pass through each intermediate length 
A < A, taking time T^^ (0, +) (note the subscript A to here!) and then reach A from there, 
i.e. 

Tn, (0, +) = Tn, (0, +) + Tn, (A, +) , (44) 

showing that the MFPT for a MT starting in the g an arbitrary length can 

be expressed fully in terms of MFPTs starting from the zero-length state. 

Although due to the fundamental asymmetry of the problem, there is no corresponding 
simple rule for the shrinking case, the following argument shows how we can leverage the 
results of the growing case to obtain a fairly compact representation. We first introduce the 
survival probability of the shrinking state with respect to rescues, which is simply given by 

(t|— ) = exp {—p^r). If no rescue occurs, the shrinking MT will hit zero length at the 
deterministic time r_ (A) = A/z/__, so that 



5, 



[T — T 



^nA(r|A,-) = (l-i7(r-r_(A)))|5_(r|-) + dr'a^{r' 
rr-ix) 

H{r-r_ (A)) { dr' a (r'|-) Sn^ (r - r'|A - z/_r', +) + (r_ (A) |-) Sn^ (r - r (A) 
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where H (x) is the standard Heavyside function. Integrating over all time yields the desired 
result 

»r_(A) 



1 r-r- ia; 

Tn^ (A, -) = — (!- e-"-^-^^^) + p_ / rfr e'^'-^Tf,^ (A - z/_r, +) + e'^'-^-^^^Tf,^ (iV) 
P- Jo 

1 ^ /o^-^'^^ra, (r|-)T^, (A - z/- r , + ) 

= (r_ (A) I-) j-l + {Tn, (A - z._r, | + 5_ (r_ (A) |-) [y^+Tn, (0,+)} 

(46) 

where (Tq^ (A — z/_r, +))(-;^ _) is the average MFPT of an MT that starts in a growing state 
after a single rescue from a shrinking state originally at length A at time 0, provided this 
happens before the shrinking state hits zero length. 

B. The splitting probabilities 

To calculate the splitting probabilities {oo\X, s) and R^ {oo\X, s), we first recall from 
the theory of Markov processes that expectation values of future events seen as functions of 



the initial time and state satisfy the backward Kolmogorov equation 16|. Specifically, any 
ruin probability R'^ (T|r, A, s), where K is some subset of fl, and C a conditioning event, 
satisfies 

^i?^ (T|r, A, .) = - (r|r. A, .) = - G,,, [i?^ (T|r, A, .')] (47) 

s' 

where the generator Gs,s' is the Hermitian conjugate of the operator (!33|) . i.e. 

C.^f"^^-^^ 7 ) (48) 

V P- -^-ax-P-J 

Since our process is time homogeneous, we can of course take the initial time to be r = 0. 
Letting our final time T — )■ oo, we see that -^^Rk (c>o|r. A, s) = -§^Rk (oo|0. A, s) = 0, so that 
splitting probabilities satisfy 

J2Gs,s' [R%{oo\X,s')] =0, (49) 

s' 

and are said to be harmonic. 
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For convenience sake we now drop the explicit mention of the final time and set 
R^{X,s) = -R^(oo|A,s) and R^{X,s) = {oo\X, s). We consider the latter splitting 
probabilities first, and write out (HSj) to obtain 



z/+^i?^ (A, +) - p+Ri (A, +) + p+R^ (A, -) = (50) 

-'^--^Ra (A, -) - P-R^ (A, -) + P-R^ (A, +) = (51) 

with the obvious boundary conditions R^ (A, +) = 1, R^ (0, — ) = 0. At this point it is 
useful to define the mean forward and backward run-lengths 

K = - (52) 

P+ 

A_ = ^ (53) 



The first equation allows us to eliminate R^ (A, 



i?f (A, -) = Ri (A, +) - K^Ri (A, +) (54) 



Insertion into the second then yields a second order equation 

^2 , d 



or equivalently 
as (|28|1 implies 



A+A_^i?f (A, +) - (A_ - A+) — i?f (A, +) = (55) 



^\r^(X,+)-±R^(X,+) = (56) 



aA2 ^ ' ' ^ dX 



We then obtain the following solutions 

= - 

(1-A+) 



^f(A,+) = 44^— H (58) 



Ra (A, -) - ,A _ (1 _ (59) 

In a fully similar manner the corresponding quantities R^ (A, s) are also readily deter- 
mined 

i^.-(A,+)^^^^^^ (60) 
Ra (A, -) - ,A _ (1 _ (61) 
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One checks that these forms satisfy the a priori requirements i?f (A, +) + (A, +) = 
(A, — ) + R^ (A, — ) = 1, which follow from the fact that the ultimate ruin of an MT 
on a finite length interval is sure. These splitting probabilities were also derived in 8[ by 
considering Laplace transforms of recurrence relations satisfied by the probability density. 

C. The conditional MFPTs 

With the splitting probabihties determined, we can directly calculate the conditional 
MFPTs by solving a time integrated form of the backward equation. Indeed, integrating 
f l47|) over the final time T , and recalling that i?^ (r|r. A, s) = 0, yields as a first step 

dT—R^{T\T,X,s) = -R^{oo\X,s) = - j dTY,Gs,s' [R^c {T\r,\ s')] . (62) 

At this point we would like to interchange the integration and the operation of G^.s', but as 
i?^ {T\t, a, s') tends to a constant for T — t- oo this is not directly possible. However, we can 
use the fact that the splitting probabilities are harmonic, i.e Gs,s' [-Rx (^1-^5 ■5')] = 0, and 
the linearity Gg^s' to obtain the identity 

- Gs,s' [RSc {T\t, a, s')] = Gs,s' [RSc (oo|A, s') - i?^ (T|r, A, s')] . (63) 

Substitution of this identity into Eq. (!62|) yields an integrable argument exactly of the form 
previously encountered in Eq. fl39l) . so that 

- / dTY,Gs,s' [i?^(T|r,A,s')] = / dTY,Gs,s' [i?^ (oo|A, s') - (T|r, A, s')] 

r POO "I 

= J2Gs,s' / rfT{/?^(oo|A,s')-i?^(T|r,A,s')} = [i?^ (oo|A, s') (A, s')] 

s' s' 

(64) 

Combining, Eqs. (!62|) and (!64|) . yields the sought after relation 

Gs,s' [R% (oo|A, s') (A, s')] = -R% (oo|A, s) , (65) 

s' 

which together with appropriate boundary conditions yields a closed form equation for the 
conditional MFPTs (A, s'). 

We now apply Eq. (!65|) to our problem, starting with the case of exiting at A we have 



15 



K-^R^A (A, +) Tf (A, +) - (A, +) Tf (A, +) + (A, -) Tf (A, -) = -i-i?^ (A, +) 



(66) 

-A_|-i?f (A, -) T| (A, -) - Ri (A, -) T| (A, -) + i?^ (A, +) (A, +) = --i?^ (A, -) 
aA p- 

(67) 

with boundary conditions i?f (A, +) Tf (A, +) = i?f (0, -) Tf (0, -) = 0. Eliminating 
R^ (A, -) Tf (A, -) and introducing the shorthand 6^ (A, +) = R^ (A, +) Tf (A, +) we find 
the following inhomogeneous second order equation 



a,.eS(A.+)-g3^eS(A.+) = 

- T—Wx^"" +' - tV (—^5 (A, +) + -fij (A, -)) s (A, +) (68) 
Using Eq. (jSlj) we can eliminate R^ (A, — ) in the inhomogeneous term, so that 

A'i (A, +) = - ^) (A, +) - f- + -) (A, +) (69) 

\X-P- X+p+J dX A+A_ Vp+ P-J 

The boundary equations for the resulting equation are 

(A, +) = (70) 

A+|-e^ (0, +) - (0, +) = --i?^ (0, +) (71) 

dX p+ 
The equation for the conditional MFPT for exiting at is the formally the same as fl68l) . 
together with boundary conditions (170|) and (17T|) . but with 6^ (A, +) and R^ (A, +) substi- 
tuted for 0;^ (A, +) and R^ (A, +) respectively. These inhomogeneous linear second order 
equations are readily solved, and we present the resultant rather unwieldy expressions in the 
appendix. As an explicit check we consider the probability of a growing MT to shrink back 
to the origin, in the absence of the boundary 

hm (A, +) = hm (A, +) ^ T (A, +) = 1 + A— (72) 

A-s>oo A-s>oo (P+^- — P-^+) 

Upon redimensionalizing, this expression is identical to the one derived earlier by Bicout 
from the full time and space-dependent survival probability [r?]. As an aside, we note that 
the result T(0, +) = 1 shows that the time-scale we have adopted is indeed that of the origin 
return time of an unconstrained MT, as already stated in Section IIIBI 
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Parameter 


Yeast 


Plant 


Growth speed ( /im/min) 


2.4 


4.8 


Shrinking speed (/im/min) 


9.6 


9.6 


Catastrophe rate (/min) 


0.3 


0.28 


Rescue rate (/min) 




0.42 


Nucleation rate (/min) 


0.15 


0.15t 


Mean length (no boundary) ( /im) 


8 


68.5 


Expected hfetime (no boundary) (min) 


4.17 


21.43 



TABLE II: Table of MT dynamical parameters for yeast and interphase plant cells, f: In the absence 
of available data we take this number equal to that of yeast. 

D. Application to biological data 



In order to get a feel what the results derived above mean in real-world terms, we apply 
them to two sets of fairly well characterized kinetic parameters for MTs, one derived from 
observations _on fission yeast [l^, and one on interphase Tobacco Bright Yellow-2 plant 



culture cells 19(|. These data sets are summarized in Table HTl 

We now confront these two types of MTs with boundaries located at 5/im, smaller than 
both mean lengths in the absence of boundaries and comparable to half the length of a fission 
yeast cell, 20/im double the mean length for the yeast MT and still significantly smaller than 
that of the Tobacco BY-2 MTs, and lOO/im of the order of typical lengthscale of a Tobacco 
BY-2 cell. We first consider the splitting probabilities (/, +) and {I, +), which we plot 
in Figure [31 We see that for the smallest boundary distance L = bfim both for yeast and 
plant MTs the probability to reach the boundary from zero length is already appreciable, 
and increases roughly linearly with starting length, consistent with it being dominated by 
uninterrupted growth. As we increase the boundary distance, the probabilties depend more 
strongly non-linearly on the starting length. This is most striking for the case of yeast at 
L = 100/im, where R^ (/, +) is essentially until the starting length is within the natural 
length I ^ 8fim from the boundary. 

Next, we turn to the conditional MFPTs (/, +) and (/, +). Here, we first need to 
take a little care, as for yeast the rescue probability vanishes (r_ = 0), so that the backward 
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Yeast 



Tobacco BY-2 





L= 100|jm 



20 40 60 





FIG. 3: Splitting probabilities Ra{1,+) (red line) and (blue line) as a function of the 

initial length / for different values of the distance L. 



runlength /_ is ill-defined. One can of course go through the procedure in Section UlI CI again 
setting p_ = at the outset in Eq. fH5]) . However, in this case the conditional MFPTs are 
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L 


Yeast 


Plant 




8.89 


3.29 


20/im 


119.1 


17.23 


lOOfim 


2.9 X 10^ 


266.55 



TABLE III: Table of the MFPT Tq^(0, +) (in min) as a function of the distance L. 



also readily determined from first principles 

(L-l) 



0:r«(/,+) 



: (/, +) 



jl+V+t) 



+ t 



L-l 

1 -e~— 



(i + ¥)) 



L-l 



(73) 



(74) 



where f[73|) follows because a non-rescuable MT can only reach the boundary without first 
shrinking away by growing towards it deterministically, and ( fMI) is obtained by averaging 
(i) the time to experience a catastrophe before reaching the boundary plus (ii) the time 
to shrink to zero length from that moment on over the ensemble of histories that do not 
reach the boundary. These two approaches indeed give the same results, which serves as 
another independent check on the general formalism. Figure H] shows the resulting passage 
times. A perhaps at first sight puzzling feature of these results is the decrease of (/, +) 
for increasing starting length I, which is evident for the yeast case. This, however, is a 
direct consequence of the conditioning on shrinking back without reaching the boundary: 
If the starting length is within the forward runlength /+ from the boundary, a conditioned 
MT must rapidly undergo a catastrophe after which it deterministically shrinks back to 
zero length. The conditioned return time (the second term on the far right hand side of 
Eq. fTMl) ) is therefore a strongly non-linearly decreasing function of the distance to the 
boundary, whereas the time to deterministically shrink back from the starting length 
only increases linearly with length. 

Finally, in Table UTTl we give the MFPT T^^ (0, +) for reaching the barrier starting from 
zero length in the growing state as calculated from Eq. f H3l) . 

We see that for a yeast MT the largest boundary distance L = lOO/im is effectively 
unbridgeable, and that even the plant MT needs on average ~ 10 times its natural lifetime 
t ~ 20 min to first reach the boundary, although it is only ^1.5 times its natural length of 
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Yeast 



Tobacco BY-2 




FIG. 4: Conditional MFPTs Tf (/, +) (red line) and (/, +) (blue line) as a function of the initial 
length I for different values of the distance L. 

I ~ 65/im. 



IV. THE GOPALAKRISHNAN-GOVINDAN SEARCH- AND-CAPTURE MODEL 



The first-passage-time model Gopalakrishnan and Govindan recently introduced ( 10| . 



hereafter referred to as GG) considers the problem of the 'capture' of a chromosome by a 
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MT 'searching' for it. It has the following ingredients: The MT is nucleated at a rate r„ 
from a centrosome in an arbitrary direction into a cone with solid angle opening of Afl. 
The centrosome is located at a distance d from the chromosome, which has a cross-sectional 
area a and therefore subtends a solid angle AQ^ = a/d'^ as seen from the centrosome. The 
probability of being nucleated into a direction in which the target can possibly be hit is 
therefore given by pc = AQc/AQ. When the MT is nucleated outside of the Harget cone\ it 
can potentially grow until it hits a cell boundary located at a distance we will call D from 
the centrosome. At this boundary, the MT is initially stalled, but experiences an increased 
catastrophe rate r?, > r^. 

We will now revisit this model, using the formalism derived in the previous sections. The 
state space of this model is conveniently represented hj Q = N U A^U B U A^U C. Here, as 
before, N is the nucleation state. At are the active states with lengths in the interval [0, D] 
in the directions that do not 'see' the target, B the state of being at the cell boundary, A^ 
the active states with lengths in the interval [0, d] and directions within the target cone, and, 
finally, C the state of being on the target chromosome. We non-dimensionalize using the 
same prescription as in Section Hi C\ denoting the additional parameters needed hj A = D/l 
and 6 = d/l. Using the results of Section HlI Al we can immediately write down an expression 
for the search-time starting from the nucleation state 

Tn/c (AT) = 1 + (1 - p,) {Ta,ub ((0, +)J + Tn/c (N) } 

Pn 

+ p, {Rl ((0, +)^J TZ ((0, +)^J + Rl ((0, +)^J {Tl ((0, +)^J + Tn/c (N)}} (75) 

The logic of this equation is simple. Starting from the nucleation state the MT (on average) 
waits 1/pn before being nucleation. With probability (1 — Pc) the nucleation will be in a 
direction that can not hit the target. In that case the MT will spend the origin-return 
time Ta^ub ((05+)yiJ this part of state space before shrinking back to zero length and 
starting again from the nucleation state. With probability Pc the initial nucleation is inside 
the target cone. In that case the MT either hits the target, without first shrinking back to 
zero-length, with probability R^^ ((0, taking time T^^ ((0, with probability 

((0) +)yic) ((0' shrinking back to zero-length before hitting the target, 

taking time ((0, +)^ ) and then trying again from the nucleation state. This process is 
illustrated in Figure O 
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A 5 




FIG. 5: Schematic illustration of the structure of the search process in the Gopalakrishnan- 
Govindan search-and-capture model. From the nucleation state the MT must either perform (with 
probability 1 — pc) a fruitless search (gray arrows) in the directions where it can interact with the 
cell boundary ^f,, or (with probability pc) it is nucleated in the right direction and then either 
directly traverses Ac to reach the chromosome C or shrinks back to zero length without reaching 
the target and must try again. 

Extracting Tq^/c (N) from the relation fl75]) yields 

pMl ((0, +)^J To/c (iV) = — + (1 - pc) Ta,ub ((0, +)J 

P+ 

+ {R% ((0, +)^J Tl ((0, +)^J + Rl ((0, +)^J Tl ((0, +)^J } (76) 

As we show in Appendix [Bl this is, apart from the changed notation, precisely the result 
derived by GG (their Eq. (34) in Appendix A) from an explicit sum-over-histories argument. 
We also note the calculation of the fruitless search time, Ta^ub ((0, +)aJ) directions 
not containing the target, is also readily simplified using the methods presented here 

Ta,ub ((0, +)aJ = Rl ((0, +)^J TZ ((0, +)^J 

+ K ((0, +)aJ [K ((0' +)aJ Ta,ub ((A, -)^J | , (77) 

where in turn 

Ta,ub ((A, -)J = Rl ((A, -)J Tl ((A, -)^J (78) 
+ Rl ((A, -)^J \tI ((A, -)^J + 1 + Ta,.b ((A, -)^J } . 

The latter expression allows Ta^ub ((A, — )^J to be expressed solely of splitting probabilities 
and conditional MFPTs. While GG use an ingenious symmetry argument interpreting a 
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shrinking MT as a growing 'anti'-MT to calculate these latter quantities, we point out 
that they can also be obtained in a straightforward manner from the differential equations 
presented in Sections llll Bl and IIII C[ 

V. CONCLUSIONS AND OUTLOOK 

Our aim was to present a structured approach to the problem of MTs interacting with 
boundaries. To this end we relied exclusively on 'backward' techniques, focussing on survival 
probabilities and their associated boundary value problems. The upshot of this approach is 
that it allows one to decompose a complex MFPT problem a priori into closed form self- 
consistency problem involving a small set of relevant splitting probabilities and conditional 
MFPTs that readily follow from a proper disjoint decomposition of the state space. The util- 
ity of this approach is illustrated by its application to the Gopalakrishnan-Govindan model, 
where the key decomposition of the search time in terms of the time spent fruitlessly search- 
ing in the wrong directions, waiting in the nucleation state and finally reaching the target is 
the starting point of the calculation, rather than, as in {lo|, the result of collecting the results 
of intermediate steps in the calculation. We hope that the technique presented will serve 
as a convenient starting point for future applications to current problems in microtubule 
cytoskeleton organization. 
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Appendix A: Explicit solutions of the conditional MFPTs 



In order to give a fairly concise expression for the conditional MFPTs we need to introduce 
a number of convenient coefficients 



a = A+ - 1 (Al) 

7 = ^-- (A3) 

P- P+ 

5 = A+{/3(l-a)+7} (A4) 
£ = 1 + a - /3p+ + a^/3p+ — 7p+ - a^p+ 
= A+ - pj (A5) 

With these definitions we find 



(A, +) = /3A+ + /3A+ (A6) 
^|-e' (2/3 + 7) (A-A) + <(A,+) (^i- (e^ - l) + 5 - a (2/3 + 7) | 

r^(A+) = -^_^^ £ + ap+(^ + 7)A + Q;/3p+A + (/3 + 7)p+(A-A)e^ . 

+ -i?^ (A, +) (£ + ap+ (/3 + 7) A + a^p+A) 

In some cases we can use the known relations 



to simplify even further. An example is the observation that 

A- \p+ P-J \K J P+ A_p_ p+ 
Taking the limit A — >■ 00 then yields, as in this limit (A, +) = 1, 

T^(0,+) = 1 (All) 

as claimed in the main text. 
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Appendix B: Formal correspondence with the Gopalakrishnan- Govindan model 



Here we provide the translation between the results of GG and our own, by noting the 
following correspondences. 



GG 


This work 


Description 


^{d,T) 


{T\ (0, +)J 


Conditional waiting time distribution reach- 
ing target without shrinking back to zero 


Qd (T) 




Conditional waiting time distribution shrink- 
ing back to zero without reaching target 


Qr{T) 


<(r|(o,+)J 


Conditional waiting time distribution shrink- 
ing back to zero without reaching boundary 




'JA.uB {T\ (0,+)^^,3T<r: = A) 


Waiting time distribution return to length 
zero after reaching the boundary at least 
once 



(Bl) 



From these correspondences we derive the identities 

poo 

3 {d, 0) = / dTal (T| (0, +)^J = Rl ((0, +)^J (B2) 

{d, 0) = - / dTTal {T\ (0, +)J = -Rl ((0, +)J ((0, +)J (B3) 
Jo 

Q'{d,0) = - / dTTal {T\ {0,+)J = {{0,+)JtZ ((0,+)^J (B4) 

POO 

Q'{R,0) = - dTTal{T\{0,+)J = -Rl{{0,+)jTg{{0,+)^^) (B5) 
Jo 

poo 

^' (0) = - / dTTaA,uB (T| (0, +)^^ , 3t < T : A, = A) 
Jo 

= -< ((A, -)J ((A, + 1 + T^,^, ((A, -)J } (B6) 
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T, 



The timescales GG introduce are therefore 

S' {d, 0) + Q' {d, 0) _ < ((0, +)^J Tl ((0, +)^J + Rl ((0, +)^J T£ ((0, +)^J 
^{d,0) i?l((0,+)J 

(B7) 

Q' iR,0) + ^' (0) 



R 



${d,0) 

Rl ((0. +)aJ Tj ((0, +)^J + < ((A, -)J {tI ((A, -)^J + i- + T^,uB ((A, -)^J 

^1 ((0, +)aJ 

^1((0,+)^J 



(B8) 



so that, indeed, their expression 



(T) = T, + ^^^T« + -n (BIO) 



fully coincides with Eq. ( !76|) . 
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